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ABSTRACT 

Using the near infrared fluxes of local galaxies derived from Cosmic 
Background Explorer (COBE)/Diffuse Infrared Background Experiment 
(DIRBE)lll J(1.25 /im) K (2.2 /xm) & L (3.5 /im) band maps and published 
Cepheid distances, we construct Tully-Fisher diagrams for the nearby galaxies. 
The measured dispersions in these luminosity-linewidth diagrams are remarkably 
small: aj = 0.09 magnitudes, ax = 0.13 magnitudes, and ai = 0.20 magnitudes. 
These dispersions include contributions from both the intrinsic Tully-Fisher 
relation scatter and the errors in estimated galaxy distances, fluxes, inclination 
angles, extinction corrections, and circular speeds. For the J and K bands, 
Monte Carlo simulations give a 95% confidence interval upper limit on the true 
scatter in the Tully-Fisher diagram of aj < 0.35 and ax < 0.45. 

We determine Milky Way's luminosity and place it in the Tully-Fisher 
diagram by fitting a bar plus exponential disk model of the Milky Way to 
the all-sky DIRBE maps. For "standard" values of its size and circular speed 
(Sun-Galactic center distance Ro = 8.5kpc and 6o = 220km/s), the Milky Way 
lies within 1.5cr of the TF relations. 

We can use the TF relation and the Cepheid distances to nearby bright 
galaxies to constrain _Ro and Qq- — log (-Ro/8.5kpc) + 1.63 log (0o/22Okm/s) = 
0.08 ± 0.03. Alternatively, we can fix the parameters of the Galaxy to their 
standard values, ignore the Cepheid zero-point, and use the Tully-Fisher relation 
to determine the Hubble Constant directly: i^o = 66 ± 12 km/s/Mpc. 

We have also tested the Tully-Fisher relation at longer wavelengths, where 
the emission is dominated by dust. We find no evidence for a Tully Fisher 
relation at wavelengths beyond lO/xm. The tight correlation seen in L band 
suggests that stellar emission dominates over the 3.3 /im PAH emission. 

Subject headings: Galaxy: fundamental parameters Galaxy: general 
galaxies: distances and redshifts 



^COBE data sets were developed by the NASA Goddard Space Flight Center under the guidance of the 
COBE Science Working Group and were provided by the NSSDC. 
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1. Introduction 

The Tully-Fisher (TF) relation between the luminosity and linewidth of spiral galaxies 
(Tully & Fisher 1977) has been used extensively as a distance indicator and to map large 
scale flows of galaxies (cf., Strauss & Willick 1995 and Jacoby et al. 1992 for reviews). Its 
usefulness as a distance indicator is limited by the intrinsic scatter of the relation. This 
scatter is lowest in redder bands where dust extinction is low: H band (1.65//m) (Aaronson, 
Huchra & Mould 1979, Aaronson, Mould & Huchra 1980, Aaronson et al. 1989, Preedman 
1990, Pierce & Tully 1992) and I band (0.90/ini) (Bernstein et al. 1994). In this paper, we 
extend the TF relation to longer wavelengths. The DIRBE experiment, with its excellent 
calibration and large beam width, is ideal for measuring the total flux of galaxies in the 
local galaxies. We describe our galaxies data set and present the results of our analysis in 
section 2. 

The Milky Way has often been deemed unsuitable for zero point calibration of the 
TF relationship, mainly because of difficulties in estimating its total luminosity. Such 
difficulties can be overcome at infrared wavelengths, where dust absorption is small. In 
section 3, we use a three-dimensional model of the Milky Way based on the DIRBE J, K, 
and L band maps to obtain a measurement of the Galaxy's luminosity. This luminosity 
can be used to place the Milky Way on the TF diagram with nearby galaxies, to constrain 
Galactic parameters, and to obtain an independent calibration of the Hubble constant. 

2. Nearby Galaxies 

For this study, we chose a sample of nearby bright spiral galaxies with sizes not too 
much smaller than the beam, measured Cepheid distances, no bright stars in a 1° field 
centered on the galaxy, and with inclination greater than 45°. We also required that the 
measured DIRBE flux from the galaxy exceed the week-to-week fluctuations in the DIRBE 
data. The galaxies used are M31, M33, M81, NGC 300, and NGC 2403. A sixth galaxy, 
NGC 247 shows a marginal detection and meets our selection criterion in only 2 bands (J 
and L). We will do our analysis both with and without NGC 247. Barring M31 and M33, 
these galaxies are smaller than the 0.7° x 0.7° DIRBE beam. M33 is slightly larger than the 
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beam, and M31 and Milky Way are extended objects compared to the beam. 

The measured fluxes were corrected for extinction using the prescription of TuUy & 
Fouque (1985), apphed to our sample using inclination angles from Pierce & TuUy (1992) 
and the standard extinction law from Mathis (1990). These corrections were never large, 
attaining maximum values of 0.18, 0.07, and 0.03 magnitudes in J, K, and L for M31. 

The distances and distance errors were taken from published literature as follows: M31 
(Freedman 1990), M33 (Freedman, Wilson & Madore 1991), M81 (Freedman et al. 1994, 
Hughes et al. 1994), NGC247 (Catanzarite, Freedman & Madore 1996). NGC300 (Freedman 
1990) and NGC2403 (Freedman 1990). Line widths for these galaxies are taken from Pierce 
& TuUy (1992) and Freedman (1990). 

2.1. Flux Extraction 

DIRBE was one of the three experiments on board the COBE mission. This experiment 
observed the entire sky in ten wavebands (from 1.25 fim to 240 fim) over a nine month 
period. Each point on the sky was observed many times, at a range of solar elongation 
angles. This gives a variable viewing geometry from within the interplanetary dust cloud, 
which is a major source of foreground emission at these wavelengths. A detailed description 
of the COBE mission can be found in Boggess et al. (1992). Details of DIRBE data 
processing, calibrations and photometry are given in the DIRBE Explanatory Supplement 
(Hauser et al. 1995). 

The fluxes of the galaxies were derived from the weekly sky maps provided by the 
DIRBE team. Intensity in each pixel of the weekly map, which is typically half a beam 
across, is a robust average of all observations in a week pointing at the region of the sky in 
that pixel. Thus the center of the beam for any observation may be up to half a pixel away 
from the center of the pixel, and does not in general lie on the center of the galaxy which 
lies in that pixel. Since the beam shape is not flat this can influence the flux measured for 
that galaxy. This error is estimated to be a few to nine percent in the DIRBE Explanatory 
Supplement (section 5.6.6). 

For each weekly map we flt a flrst order polynomial to the sky brightness in a small 
annulus around the pixel that contains the galaxy. The sky background was then subtracted 
from the maps. A point source can influence the flux levels of neighboring pixels due to the 
way the weekly maps were made, so pixels adjacent to galaxies were excluded from the sky 
estimate. All the weekly maps which had observations of that part of the sky were then 
averaged (after rejecting outliers) to give the the average flux of the galaxy. The week to 
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week variation of the flux in a given pixel gives an estimate of the noise in the maps. The 
flux per pixel thus obtained in MJy/Sr is multiplied by the beam size (given in the DIRBE 
Explanatory Supplement) to derive the flux of the sources. 

The flux in a given pixel has contributions from the galaxy, unresolved Galactic stars 
and zodiacal hght. The zodiacal hght varies from week to week as one sees the same part 
of the sky through a different path length of the zodiacal dust cloud near the Earth. This 
makes it necessary to remove the background separately for each week. We assume that 
the starlight is smoothly varying in space and can be subtracted as part of the background. 
This is true except in the rare case of a very bright star close to the galaxy. Wc examined 
the Palomar Observatory Sky Survey plates for such bright stars near candidate galaxies, 
and rejected NGC 925 and NGC 4571 from our sample because of nearby bright stars. The 
week-to-week variation in the flux of the source gives the estimate of errors due to imperfect 
zodiacal light subtraction, and also due to the beam being centered at different positions 
with respect to the galaxy. We have used the week-to-week flux variation as a conservative 
estimate of the errors in the flux measurement. 

M33 is slightly bigger than the beam size, and the single pixel method described above 
underestimates its flux. We simulated observations of M33 assuming an exponential disk 
with a scale length of 6' (Regan & Vogel 1994) and the beam profile of the DIRBE beam. In 
these simulations, we placed the beam randomly within the pixel closest to the position of 
M33 center. The same experiment was repeated with a point source. The flux measured for 
the M33 simulation was compared to the flux measured for a point source and a correction 
factor of -0.55 magnitude was added to the M33 magnitude. M31 appears extended in the 
DIRBE maps, so its flux was measured by summing all the pixels containing the galaxy 
and normalizing this sum by the ratio of the pixel area to the beam area. 

2.2. Mctximum Likelihood Analysis 

Wc determined the TF parameters and errors using a maximum likelihood analysis. 
This approach is preferable to a fit to the data as the maximum likelihood analysis 
includes the selection bias due to our magnitude limit (Willick 1994). In a magnitude hmited 
sample, selection bias implies that galaxies near the luminosity cut-off are systematically 
brighter (Willick 1994). NGC 247, which falls just below our cut-off, behaves just as 
expected given this effect: it is brighter than the best fit TF relation. We perform our 
analysis both with and without NGC 247. 
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The maximum likelihood function for our sample is, 

P{rh\a,b,aTF) = Yl 
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Here the observed quantities are the galaxies' magnitudes m^, Cepheid-based distance 
moduli Di, the logarithms of their line widths rji, and the magnitude limits mi^i (i.e., the 
amphtude of the week-to-week fluctuations measured at the galaxy location) . The model 
parameters a, b, and axF are the slope, zero-point, and [total] scatter about the TuUy-Fisher 
relation. Note that axF includes contributions from both the intrinsic TuUy-Fisher relation 
scatter and the errors in estimated galaxy distances, fluxes, inclination angles, extinction 
corrections, and circular speeds. The product is evaluated over all galaxies in the sample, 
i.e., all galaxies for which mj < m^^j. 

The maximum likehhood flts to the J, K and L band TuUy Fisher relations (Figures 1, 
2 & 3) are: 

r = -8.13l^;l^log(A\//400km/s) - 22.00l^;i^ 

= -8.59lg:^^log(Ay/400km/s) - 23.0llg:Ji 

= -9.011°:^^ log(Al//400km/s) - 22.991°;}^ 

where J'^, K'^ and are dust-corrected magnitudes. The linewidth-luminosity, or TuUy 
Fisher, relation in these bands is very tight, showing a scatter of (Ttf{J) — 0.09 magnitudes 
in J band, (Ttf{K) — 0.13 magnitudes in K band and (7tf{L) — 0.20 magnitudes in L 
band, flts yield smaller (but less behevable) values for a. The minimization of the 
likelihood function was done with the Numerical Recipes (Press et al. 1992) AMOEBA 
program. Investigation of the likelihood surface showed that it was well-behaved near its 
minimum. A small scatter in the TF relation of 5 galaxies implies a small intrinsic scatter 
in the TF relation for all galaxies. We confirm this for our sample with a Monte-Carlo 
simulation with 1000 realizations of the local galaxies sample. We generated 5 galaxies with 
the reported 77,, and our best fit slopes. We varied the value of (Ttf in our simulations. 
For a true aTF{J) — 0.35, we recovered a scatter as small as observed in less than 5% of 
the simulations. For axpiJ) — 0.55, we recovered such a small scatter in only 1% of the 
simulations. The quoted errors on the TF slope and zero-point are the 95% confldence 
intervals from the Monte-Carlo simulations. 

If we include NGC 247 in the sample, then the maximum likelihood flts are 

r = -7.781og(Al^/400km/s) - 22.07 , 
= -8.901og(Ay/400km/s) - 22.99 , 
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8.761og(Ay/400km/s) - 23.00 , 



and the corresponding Tully Fisher relation scatter is aj — 0.15, ax = 0.14, and = 0.21 
magnitudes. For a true arriK) — 0.45, we found such small scatter in less than 5% of 
simulations. For (Ttf{K) — 0.7, we found such smaller scatter in only 1% of simulations. 

The dust emission at longer wavelengths 60-240 /xm does not show a luminosity- 
hnewidth relation. From this and from the tight L band Tully Fisher correlation, we 
conclude that the non-stellar or dust contribution to the L band is not large. This is 
consistent with recent estimates of 8-16 % contribution of the 3.3 //m PAH feature to the L 
band (Bernard et al. 1994) 

While giant stars are thought to contribute most of the light in near IR bands, 
supergiants can make an important contribution in the spiral arms (Rhoads 1996). This 
contribution does not appear to be a source of significant scatter in the TF relation. 



We cannot directly use the total flux from the Milky Way (MW) as regions at different 
distances contribute to that flux. We estimated the luminosity density of the MW from the 
best fitting model to the DIRBE data. The MW was modelled as a sum of an exponential 
disk and a bar (Spergel, Malhotra & Blitz 1996, Paper I). Extinction corrections were based 
on a 3-dimensional dust model. Varying modelling parameters and extinctions changes 
the MW flux by less than 10a lOcircular speed Oq = 220 ± 10 km/s and a distance to the 
Galactic center Rq = 8.5 ± 0.5 kpc (Gunn, Knapp & Tremaine 1979). Using these values, 
we obtain absolute magnitudes J = —23.05, K = —24.06, and L = —23.88. More recent 
determinations of the Sun-Galactic center distance using water masers suggest a smaller 
distance: 7.1 ± 1.5 kpc (Reid et al. 1988) and 8.1 ± 1.1 kpc (Gwinn, Moran & Reid 1992). 
Combining our uncertainties and assuming a TF slope of 9 implies an intrinsic uncertainty 
in the Galaxy's position in the TF diagram of 0.23 magnitudes. 

If we include Milky Way in the sample, then the maximum likelihood fits are: 



3. The Milky Way 



r = -8.091og(Ay/400km/s) - 22.15 



-9.231og(A\//400km/s) - 23.03 
8.911og(Al^/400km/s) - 23.04 



For uniformity, the velocity width of the MW was scaled from the observed velocity width 
of M31: /W{MW) = (220/250) x A\/(M31). The Tully Fisher relation in these bands 
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remains very tight with aj — 0.18, gk — 0.16, and gl — 0.21. We repeated the Monte 
Carlo simulations as described in Section 2. For a true axriK) = 0.45, we found a scatter 
as small as observed in less than 5% of simulations. For axriK) = 0.57, we found such a 
small scatter in only 1% of simulations. With the standard parameters, the Milky Way is 
~ 1.5(7 too luminous in J, ~ la too luminous in K, and ~ 0.3a too luminous in L. This 
basic agreement shows that the Cepheid distance scale is consistent at the 10% level with 
Gunn et al. (1979) distance determination to the Galactic center. 

Alternatively, wc can use the TF relation obtained for the 5 galaxy sample, together 
with our model for the Galactic luminosity to constrain Galactic parameters: 



-Slog (-4^ ] +8.13 log f — =0.42 ±0.16 



^8.5kpcy ■ ''V220km/s^ 

-5 log I ^° I +8.59 log I — =0.39 + 0.19 
^\8.5kpcy ^V220km/sy 

-5 log I — ^ 1+9.01 log I — 1=0.20 + 0.24 
^\8.5kpcJ ^V220km/sy 

where the 3 equations correspond to fitting the J, K and L band relations respectively. 
The uncertainties include the width of the TF relation, the quoted uncertainties in the 
Galactic model (10%) and a conservative estimate of systematic errors in point source 
flux measurements (9%). This comparison of the Galaxy with the nearby galaxies sample 
suggests that its rotation speed is unlikely to be as small as 200 km/s. 

Finally, we can also use the Galaxy (if we are willing to assume values for Rq and 
©o) to obtain a value of the Hubble parameter that is independent of the zero-point of the 
Cepheid distance scale (Wright 1994). If we use the Galaxy as the zero-point of the K band 
TuUy Fisher relation, we find that 

Rn \ / ©0 \ 



= -8.59[log(Ay) - 2.5] - 22.52 + 0.29 - 5 log — + 8.59 log 



8.5kpcy ■ ''V220km/sy 

where the zero-point error is the sum in quadrature of the uncertainties in Galactic 
parameters, the systematic error in the point source flux measurements, the uncertainties 
in the Galaxy model, and the deviation of a typical galaxy from the TF relation. 

In order to extrapolate to H band, where there are measured fluxes for large numbers 
of external galaxies (Aaronson et al. 1986), we use measured H band fluxes for our sample 
(Freedman 1990, Pierce & TuUy 1990) to determine the dust-corrected H'^ — K"^ and 
jc _ color-color relation for the galaxies in this sample. A flt to the data yields 
{H" - K") = 3.27 - 2.50 x {J" - K") with a scatter of 0.18 magnitudes. This suggests 
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{H^ - K") = 0.75 ± 0.18 for the Milky Way and therefore, H" = -21.77 ± 0.34 for a galaxy 
with log(A?,') = 2.5. If we use this estimate to set a zero-point for the H band TF relation 
and combine it with the Aaronson et al. (1986) TF study of the local velocity field (their 
Table 6 and Equation 4), this yields i^o = 66 ± 12 km/s/Mpc. This value is lower than that 
obtained by using Cepheid distances and the Local Group to normalize the TF relation 
(Aaronson et al. 1986, Freedman 1990) as the Galaxy is somewhat over-luminous relative 
to the rest of the nearby galaxies sample (for the standard Galactic parameters and the 
standard Cepheid distances). If we assume that Rq — 7.1 kpc (the smaller maser-based 
measurement [Reid et al. 1988]) and hold ©o constant, then the Hubble constant estimate 
increases to 79 km/s/Mpc. If we assume that Rq — 1.1 kpc and hold the angular rotation 
rate, 0o/-Ro constant, then the Hubble constant estimate decreases to 58 km/s/Mpc. 
Besides the uncertainty in Galactic parameters, the dominant source of error in these 
estimates is the extrapolation from J and K band fluxes to H band. Thus, a K band survey 
of spirals in the Goma cluster would yield a more definitive value of the Hubble constant. 

4. Conclusions 

We find that TF relation extends to longer wavelengths (2.2 ^m and 3.5 /xm) than 
previously explored. The extinction corrections (cf. Mathis 1990) at 2.2 /im and 3.5;um 
are about half and one-third as large as for the H-band at 1.65/im, so these bands may be 
usefully exploited for distance estimation. With the advent of imaging IR instruments and 
two major near-IR sky surveys (2MASS and DENIS), there is also potential to estimate the 
distances and hence the peculiar velocity flows for many more galaxies ( ~ 10^) in a greater 
part of the sky and nearer to the plane of the Milky Way. 

Using Gepheid distances for galaxies in the Local Group (not including the Milky Way) 
and assuming standard Galactic parameters, we flnd that the Galaxy obeys the TF relation 
for nearby galaxies with Gepheid distances. This consistency is an independent check of the 
distance scale. 
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Fig. 1. — The J band (1.25/im) TuUy-Fisher diagram. The error bars for the local galaxies 
M31, M81, NGC 2403, M33, NGC 300, and NGC 247 represent the la uncertainties 
in the dust-corrected absolute J band magnitude (J'^) due to photometric errors and 
reported uncertainties in Cepheid distances, added in quadrature. The Milky Way (MW) 
absolute magnitude assumes a Galactocenteric distance of 8.5 ± 0.5 kpc and circular speed 
Go = 220 ± 10 km/s; its error bar represents the 0.23 mag uncertainty derived in section 3. 
The solid line is the fit to the 5 galaxy sample, excluding NGC 247 (which falls just below 
our magnitude limit) and the Milky Way. The dashed line is the fit to the full 7 galaxy 
sample. 

Fig. 2. — As figure 1, but for K band (2.2/im). 
Fig. 3. — As figure 1, but for L band (3.5/im). 



